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[Bi(Zni/2Tii/2)03]^. [PbTi03]i_^ (BZT-PT) solid solution. We find that hybridization be- 
tween Zn 4p and O 2p orbitals allows the formation of short, covalent Zn-0 bonds, enabling 



'^ ' favorable coupling between A-site and B-site displacements. This leads to large polarization, 

O ■ strong tetragonality and an elevated ferroelectric to paraelectric phase transition temperature. 

Inhomogeneities in local structure near the 90° domain boundaries can be deduced from the 

r> I asymetric peak broadening in the neutron and x-ray diffraction spectra. These extrinsic effects 

S ■ 

make the ferroelectric to paraelectric phase transition diffuse in BZT-PT solid solutions. 



Due to their importance in ferroelectric and piezoelectric applications and their fas- 
cinating physical properties, perovskite solid solutions have been an area of intense re- 
search [l|, I2, S y, |5| . Recently, it was shown that very large c/a ratios (1.11) and high Tc 
(^1000 K) can be obtained in the [Bi (Zui/aTii/a) Og]^. [PbTi03]i_^. (BZT-PT) solid solu- 
tion 6]. This makes BZT-PT a promising candidate material for use in next-generation 
piezoelectric and nonvolatile memory applications. In particular, the rate of increase in 
both properties as BZT replaces PT is the largest known for any PT-based perovskite solid 
solution. Experimental investigations of dielectric properties in the BZT-PT system also 
show that the ferroelectric to paraelectric phase transition becomes increasingly diffuse with 
greater BZT concentration. This is unusual, as increased tetragonality and stronger fer- 
roelectricity is typically accompanied by a sharpening of the dielectric constant peak at 
the Tc and loss of relaxor and/or diffuse phase transition behavior [7]. We note that the 
Bi(Mg,Ti)03 - PT (BMT-PT) solid solution, which does not exhibit such anomalous behav- 
ior [s|, differs from the BZT-PT system only by the isovalent B site cation replacement of 
Mg for Zn. In this Letter, we combine theoretical calculations with new experimental data 
to reveal the origins of anomalous behavior in the BZT-PT system. 

To determine the local, A-scale structure in the BZT-PT system, we performed density 
functional theory [9| (DFT) calculations using 2 x 2 x 2 40-atom supercells with periodic 
boundary conditions at experimental volume [6[. Technical details are the same as in pre- 
vious work |li|. We use five and seven different cation arrangements with minimal oxygen 
over- and „„der-bo„d,ng Qy Q to .tudy .^0.25 and .^0.5 compositions respectively. 
Polarization magnitudes were calculated using the Berry-phase approach J13l ]. 

Results of theoretical calculations were compared with experimental structural informa- 
tion obtained from neutron powder diffraction (NPD) of selected compositions in the (x)BZT 
- (l-x)PT system. The neutron powder diffraction data were collected under ambient condi- 
tions over the range 3-165° 29 using BT-1 32 detector neutron powder diffractometer at the 
NIST Center for Neutron Research. Full details of the experimental data and the Rietveld 
structural refinements are available as supplementary online material. Additional structural 
information was extracted from the neutron scattering data collected at NPDF spectrom- 
eter at LANSCE, Los Alamos National Lab by pair distribution function (PDF) analysis 
techniques; details of this experimental procedure have been described previously |l4| . 

A summary of structural data obtained from both DFT calculations and refined NPD 



data are shown in Table I. This compares the measured and predicted (100) cation dis- 
placements for several compositions in the (x)BZT - (l-x)PT system and reveals very good 
agreement between the theoretically and experimentally determined values. While displace- 
ments determined by refinement of the NPD data show only average A and B site behavior, 
in the DFT calculations it is possible to distinguish between the chemically distinct cations 
occupying each cation site in the tetragonal perovskite structure. A closer examination of 
the different behavior of the cations in the (x)BZT - (l-a;)PT system helps us understand 
its anomalous behavior. 

We find that Zn distortions are very large, reaching 0.5 A in the x=0.50 composition. 
The combination of the low valence of Zn and the long, ^2.6 A Zn-0 distances created 
by the large displacements eliminates the bond between Zn and one of its six O neighbors, 
transforming the Zn-Og cages are into Zn-Os complexes. The unusually large Zn distortions 
are due to the hybridization of 4s and 4p orbitals of Zn with 2p orbitals of O atoms. Local 
atom projected density of states (LDOS) (Figure la) shows that Zn 4s and 4p orbitals are 
partially filled in (x)BZT-(l-a;)PT, due to covalent bonding with oxygen 2p states. This is 
similar to the LDOS of covalently bonded Ti ions in PbTiOs (where the hybridization is 
between Ti 3d and O 2p) jl,], and is in stark contrast to the highly ionic bonding of Mg 
ions in BMT-PT J29|. In comparison with Zn, Ti displacements in the (x)BZT - (l-x)PT 
system do not change with alloying x. From a;=0.25 to a;=0.50, the calculated Ti cation 
displacement only increases from 0.32 A to 0.34 A. Even for a;=0.50 composition, the Ti- 
Oq cage is largely preserved, with average bond order RiO.25 [l^ for the long Ti-0 bond 
distances. The large Ti displacements are enabled by the Zn off-centering, due to coupling 
between B-cation displacements |l0lll6l|. 

Both our experimental and theoretical data (Table 1) show that large Bi and Pb dis- 
tortions are present in the BZT-PT solid solutions. The typical Bi distortion magnitudes 
are about 0.80 A and 0.90 A for x=0.25 and x=0.50, respectively. The magnitudes of Bi 
displacements found in this study are much larger than the 0.5 A distortion found for Bi in 
BiFeOs |17i]. Addition of BZT also induces a 10-20% increase in Pb displacement magni- 
tudes. Analysis of Bi-0i2 and Pb-0i2 cages reveals a crucial difference between the bonding 
behavior of the two A-site cations. Bi displacements create three or four very short (?^2. 1-2.2 
A) and strong (bond order ^0.6-0.8) Bi-0 bonds. This motif was also found by a recent x- 
ray absorption fine structure study of the (K,Bi)Ti03-(Na,Bi)Ti03 solid solution [isf. Bond 



valence analysis [15| shows that typically 70-90% of the Bi-0 bonding is with these three 
or four closest O atoms, while the other O atoms are located far enough from Bi to make 
their contribution to the bonding negligible. Essentially, the large Bi displacements create 
a Bi-Os or Bi-04 complex out of the Bi-0i2 cage. In contrast, although the Pb atoms also 
make three or four short Pb-0 bonds, the larger ionic size of Pb and its smaller displacement 
magnitude mean that the differences among Pb-0 bond lengths are smaller. The shortest 
three or four Pb-0 bonds typically account account for only 40-60% of the total Pb-0 bond 
order. Thus, the difference in ionic sizes of Pb^+ (1.49 A) and Bi^"*" (1.36 A) gives rise to 
different A-0 bonding motifs. 

In agreement with DFT results, we find that substitution of BZT for PT leads to a broad 
pair distribution function (PDF), indicative of a material with strong local distortions from 
overall structure. Comparison of the neutron-scattering PDF for a;=0.20 and the relaxed 
DFT PDF for x=0.25 shows excellent agreement between theory and experiment (Figure 
lb). The shoulder at 2.2 A in the experimental PDF is due to the short Bi-0 distances 
created by the large Bi displacements. In the higher-resolution DFT PDF, another shoulder 
is apparent at 2.5 A, consistent with the off-center displacement of Pb ions (l2Lll9|. 

The local structure distortions determined by our DFT calculations manifest themselves 
in the observed structural and dielectric properties of the BZT-PT solid solution. The 
increase in the c axis lattice parameter is driven by the large magnitude of Bi distortions 
due to coupling of A-site displacement energy stabilization to the strain in the material [l|, 
120, |2l|, |22|. An increase of c parameter also alleviates large direct and oxygen- mediated 
A-B repulsions created by the large Bi displacements Isl, LLO]. In contrast to other solid 
solutions, the B sites of BZT-PT are completely occupied by ferroelectrically active (Zn 
and Ti) ions that can form short, covalent bonds with oxygen. For these cations, the large 
3-site displacement necessitated by large Bi displacements and A-B displacement coupling 
are favorable. For these cations, there is no energetic penalty for making large 



3-site disp l 



off-center displacements necessitated by a combination of large Bi-cation displacements and 
oxygen bond order conservation. This stabilizes the tetragonal perovskite phase |22| . While 
addition of BZT increases the c axis parameter, the a axis parameter is virtually unchanged 
with composition. This is most likely due to the fact that Ti-0 octahedra are rather stiff, 
and there is a si gnifi cant energy cost for compressing Ti-0 bonds along the short (010) and 
(001) directions 21*. The large ionic size of Zn"*"^ also contributes to preventing decrease of 



le 



the a parameter. 

The local structural distortions (reported above) for BZT-PT also help rationalize t 
dielectric properties. Reported dielectric permittivity for the (x)BZT - (l-x)PT system 
shows that increased BZT content causes a nonlinear rise in the temperature of the dielectric 
maximum (T^^max)- The trend in Tg^ax can be understood by considering the changes of 
polarization magnitude with BZT content in conjunction with the previously established 
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2J] . Using Berry- 



relationship that Te^max is proportional to P^ for PT-based solutions 
phase calculated P values for pure PT (x=0.0), a;=0.25 and x=0.5 compositions, we see that 
the increase in T^^max correlates with the enhancement of P, with extremely high P values of 
1.03±0.05 C/m^ and 1.27±0.05 C/m^ for the x=0.25 and x=0.50 compositions respectively 
(Table I). This is due to the large cation displacements present in these compositions. 

The increase of T^^max in BZT-PT is accompanied by a greater dispersion of the phase 
transition [y]. It is unlikely that the appearance of diffuseness in the phase transition is 
due to the changes in local structure. Analysis of relaxed DFT coordinates shows only 
small (~15°) scatter in directions of cation displacements. Such small scatter indicates that 



all the cations experience nearly the same loca 
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otential energy surface; this is typical of 
2611 that do not exhibit diffusive or relaxor 



highly tetragonal PT-rich solid solutions 

dielectric behavior \il. We therefore conclude that the diffuse phase transitions observed in 

BZT-PT are not due to intrinsic properties but must be due to extrinsic effects. 

Information about the microstructure of ceramic compositions in the (a;)BZT - (l-a;)PT 
system can also be extracted from the shape of the diffraction peaks in the experimental 
NPD data. In Figure 2a, we present the the 40-60° 29 neutron diffraction data for the 



x=0.35 composition J30|. Arrows point out clearly visible asymmetric peak broadening on 
the high 26 angle side of the /-containing peak reflections ((001), (102), etc.). Although less 
marked, asymmetric broadening is also present on the low 29 angle side of /i/fc-containing 
peak reflections. It is well known that peak broadening in diffraction experiments can have 
various causes, including small crystalline particle size or an internally stressed sample. 
However, both of these possible origins would result in symmetrical peak broadening. In 
contrast, the broadening observed in the NPD presented here, and in XRD data of the 
(x)BZT - (l-x)PT system |6|, is quite strongly asymmetrical. 

We propose that this unusual broadening may be explained by a microstructural model 
based on local suppression of polarization and tetragonality near ferroelectric domain bound- 



aries in highly tetragonal compositions with large P. In non-poled polycrystalline samples 
with randomly oriented domain polarizations, large elastic strain conditions make it energet- 
ically unfavorable for domains with non-aligned polarizations to lie adjacent to one another 
across a domain boundary. However, the strain can be minimized by the creation of a nar- 
row region of locally suppressed polarization and tetragonality near the domain boundary, 
illustrated schematically in Figure 2b. In the strained near-boundary region, the a and c 
lattice parameters are distorted, with nearly cubic unit cells at the domain boundary. When 
analyzed by diffraction, a sample possessing the microstructure shown in Fig. 2b will ex- 
hibit multiple sets of peak reflections. In addition to the intrinsic reflections resulting from 
material in the center of the domains, the /-containing peak reflections of a sample with 
this structure will also exhibit a broad asymmetrical tail on the high 26 angle side of the 
fundamental peak. This broadening effect is due to reflection contributions from reduced 
c axes of the unit cells near the domain walls. The increase in the a and b lattice param- 
eters in the near-boundary region driven by conservation of unit cell volume and rotation 
of polarization [22| will give additional asymmetrical tails on the low 29 angle side of the 
fundamental /i/A;-containing peak reflections. 

To approximate the strained structure with continuously varying lattice parameters near 
the ferroelectric domain boundaries in compositions of the (x)BZT - (l-x)PT system, we 
used a two-phase model in the Rietveld refinement of the neutron diffraction data for the 
x=0.20 and x=0.35 compositions. Full details of this refinement are available online. The 
model consists of high and low c/a ratio tetragonal phases. The high c/a phase represents the 
undistorted intrinsic structure at the center of domains, while the low c/a phase represents 
an average of the unit cells in the region of reduced polarization near the domain boundaries. 
The estimated volume fractions of high and low c/a phases are 0.90 and 0.10 respectively. 
As shown in Figure 2b, the Rietveld refinement results using this model agree with the 
experimental data. The ionic displacements obtained by the Rietveld refinement for the 
high c/a phase agree with our DFT calculated values (Table I). Unsurprisingly, the low c/a 
ratio values (~0.5 A for the A-site displacement and ?^0.25 A for the B-site displacement for 
x=0.20 and a;=0.35) are quite different from DFT results, since the strain gradients present 
in the proposed microstructure were not included in the DFT calculations. 

We believe that this non-uniform local environment in the (a;)BZT - (1-2;)PT system, 
rather than any intrinsic distortions, is the cause of the increased diffuseness of the dielectric 



transition observed for compositions with greater BZT content. Local polarization varies 
between the center of each 90° domain and the domain boundary. This inhomogeneity leads 
to a distribution of local transition temperatures and a diffuse dielectric constant peak J28| . 
Additionally, the presence of uncompensated charge at the 90° domain boundaries gives rise 



to electric fields in the material [27[. It is likely that these also contribute to the smearing 
of the phase transition. Comparison of low and high cja displacements shows that as BZT 
content is increased, the differences between the domain interior and boundaries become 
more pronounced. It is this effect that leads to increased transition diffuseness. 

In conclusion, we have examined the enhanced tetragonality BZT-PT solid solution with 
a combination of experimental and first-principles methods. Favorable coupling between A- 
site and B-site cation displacements gives rise to large polarization magnitudes. Analysis of 
neutron diffraction data reveals that the high cja ratio and P magnitude lead to significant 
changes in local structure near the 90° domain boundaries. In turn, the inhomogeneity of 
the local structure smears out of the ferroelectric to parelectric phase transition. 
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TABLE I: BZT-PT structural data obtained from DFT relaxed coordinates for x=0.25 and x=0.5 
compositions and from Rietveld refinement of experimental neutron difraction data for x=0.20 and 
x=0.35 compositions. Cation displacements from center of oxygen cage in A. Cation displacement 
scatter away from (100) direction (6) in degrees. For experimentally investigated compositions, 
only average characteristics of A- and B-sites are available. Data are for the intrisic high c/a 
structure, unless noted otherwise. Low c/a phase values are for the region next to domain boundary. 
Polarization values in C/m^. 
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DFT 0.49 0.82 
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0.29 
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0.50 


DFT 0.56 0.90 


0.73 


0.50 0.34 


0.38 
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FIG. 1: (a) Local density of states for Zn and O in x=0.25 BZT-PT. A cutoff radius of 2 a.u. was 
used to perform th e projection on atomic orbitals of Zn and O. Strong hybridization is present 
between Zn 4p and O 2p orbitals. This allows formation of short covalent Zn-0 bonds and favors 
large Zn off-centering, (b) (x)BZT-(l-x)PT pair distribution functions (PDF) obtained by neutron- 
scattering for x=0.20 (solid) and computed from relaxed DFT structures for x=0.25(dashed). The 
experimental and theoretical data show excellent agreement. 
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FIG. 2: (a) Experimental and refined neutron diffraction data plotted in the 40-60° 26 range for the 
x=0.35 composition. Crosses indicate experimental data, two-phase fit shown by solid line and fit 
difference is marked at bottom. Black arrows point to asymmetric broadening of /-containing peaks, 
while the gray arrow highlights the more subtle h, k peak broadening, (b) Schematic illustration of 
locally suppressed tetragonality and polarization in the region near a 90° domain boundary. Shaded 
area indicates the strained boundary region with suppressed polarization and tetragonality. 
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